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ABSTRACT 

We use semi-analytic techniques to evaluate the burst sensitivity of designs for 
the EXIST hard X-ray survey mission. Applying these techniques to the mission 
design proposed for the Beyond Einstein program, we find that with its very 
large field-of-view and faint gamma-ray burst detection threshold, EXIST will 
detect and localize approximately two bursts per day, a large fraction of which 
may be at high redshift. We estimate that EXISTs maximum sensitivity will 
be ~ 4 times greater than that of Swiff s Burst Alert Telescope. Bursts will be 
localized to better than 40 arcsec at threshold, with a burst position as good as 
a few arcsec for strong bursts. EXISTs combination of three different detector 
systems will provide spectra from 3 keV to more than 10 MeV. Thus, EXIST will 
enable a major leap in the understanding of bursts, their evolution, environment, 
and utility as cosmological probes. 
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1. Introduction 

In its quest to find black holes throughout the universe, the Energetic X-ray Imaging 
Survey Telescope (EXIST) will detect, localize and study a large number of gamma-ray 
bursts, events thought to result from the birth of stellar-mass black holes. We present the 
methods used to calculate EXISTs capabilities as a gamma-ray burst detector; we use the 
EXIST design evaluated by the National Research Council's 'Committee on NASA's Beyond 
Einstein Program: An Architecture for Implementation' (2007; see also Grindlay 2007). The 
combination of large detector area, broad energy coverage, and wide ficld-of-view (FOV) will 
result in the detection of a substantial number of bursts with a flux distribution extending to 
fainter fluxes than that of previous missions. Thus EXIST should detect high redshift bursts, 
perhaps even bursts resulting from the death of Pop III stars. EXISTs imaging detectors 
will localize the bursts, while the combination of detectors, both imaging and non-imaging, 
will result in well-determined spectra from 3 keV to well over 10 MeV. 

In this paper we first describe the EXIST mission design (§2), emphasizing aspects 
relevant to burst detection. Then we present the sensitivity methodology (§3), which we 
apply to the individual coded mask sub-telescopes (§4). EXIST will consist of arrays of these 
detectors with overlapping FOVs, and the overall mission sensitivity results from adding 
the sensitivity of the individual sub-telescopes (§5). Imaging using counts accumulated 
over different timescales increases the sensitivity (§6). Finally, we combine these different 
calculations to evaluate EXISTs overall capabilities to study bursts (§7). 

2. Overview of the EXIST Mission 

The EXIST design analyzed here was proposed as the Black Hole Finder Probe for 
NASA's Beyond Einstein program. In this design, described in Grindlay (2007), the mis- 
sion consists of two arrays of sub-telescopes. The 19 High Energy Telescopes (HETs) will 
use a Cadmium Zinc Telluride (CZT) detector plane, while the 32 Low Energy Telescopes 
(LETs) will use a silicon detector plane. The spacecraft will be launched into low Earth orbit 
(~500 km) by either a Atlas V-551 (for an orbital inclination of i ~ 20°) or a Delta IV 4050H 
{i ~ 5°) for a 5 (minimum)-lO (goal) year mission. The Csl active shielding for the CZT 
detectors will also be instrumented to provide spectral coverage for gamma-ray bursts at 
higher energies. Table 1 provides the detector parameters relevant to this study. The space- 
craft pointing will rock ±15° perpendicular to the orbit around the zenith, resulting in nearly 
uniform sky coverage and sensitivity. 

Both the HETs and LETs will image the gamma-ray sky using the coded mask technique. 
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The detector plane 'sees' the sky through a mask with open and closed cells that is a fixed 
distance above the detector plane. Therefore a source in the FOV casts a shadow with the 
mask's pattern on the detector plane. The distribution of sources on the sky is deconvolved 
from the counts detected by the position sensitive detectors. Sources in the central part of 
the FOV, the 'fully coded' region, illuminate the full detector plane, while sources further 
out in the FOV, the 'partially coded' region, illuminate only a fraction of the detector plane. 
The dimensions of the detector plane and mask, and the detector-mask distance, determine 
the FOV, while the detector-mask distance and the dimensions of the mask cells and detector 
pixels fix the angular resolution. 

As for Swift and GLAST, EXIST will run burst detection and localization software 
onboard (the Fast Onboard Burst Alert System — FOBAS), and telemeter data to the ground 
for further analysis. In the current design, FOBAS will run both rate and image triggers 
on the datastream from both the HET and LET sub-telescopes. Rate triggers will search 
for statistically significant increases in the count rates from the sub-telescopes. The image 
triggers will form images from the counts from the individual sub-telescopes, add the images, 
and search the resulting sky image for a new, statistically significant point source. In the 
current design, images will be formed with 3, 18, 108, 648 and 1296 s accumulation times. 
When a burst is detected, EXIST will downlink the burst time and location (as well as other 
basic burst parameters) through the Tracking and Data Relay SateUite System (TDRSS) 
within ~ 10-20 s, as is done for Swift and will be done for GLAST. 

Data indicating the time, energy and pixel of every HET and LET count, as well as 
other standard science and housekeeping data, will be downlinked approximately every four 
hours through a TDRSS Ku band link, as will be done for GLAST. Ground software will 
calculate more accurate sky positions and other parameters (e.g., durations and spectra) for 
the bursts detected onboard, and will search the datastream for bursts that FOBAS did not 
detect. 

The active Csl shields behind the HET detector plane and in the lower part of the 
HET coUimators will be instrumented to provide 64 channel spectra between ~300 keV to 
~10 MeV (see Garson et al. 2006a). The current plan is that spectra accumulated every 
1 s will be downlinked. By buffering the counts from the shields, the count binning will be 
increased to every 0.1 s for the time period 500 s before to 500 s after the trigger. 
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3. Burst Detection Sensitivity 



A burst will be detected by EXIST when a statistically significant new source is found 
in either an HET or LET image of the sky. The same criterion applies to Swiff s BAT (a 
single CZT coded mask detector) and therefore the sensitivity analysis we use follows the 
methodology applied to the BAT in Band (2006). 

Formation of the image in which the burst is detected may be initiated by either a rate 
or image trigger. A rate trigger will search the count rates from the sub-telescopes for a 
statistically significant increase. An image trigger will search for new sources in images of 
the sky that will be formed continuously. The new source in an image trigger may not be 
statistically significant, but will indicate that a burst may be in progress. After either a rate 
or image trigger, FOBAS (the burst fiight software) will vary the time and energy ranges 
over which counts are accumulated to maximize the signal-to-noise ratio. An image will then 
be formed from these counts. The threshold for these initial rate or image triggers will be 
sufficiently loose to allow many triggers; the absence of a statistically significant point source 
in the final image will weed out the false positives. Note that the final image in which the 
point source is most significant may be formed from the counts in a different energy and 
time bin than that of the counts that initially triggered FOBAS. 

Regardless of the process leading to the final image, EXISTs burst sensitivity will be 
the minimum burst fiux that results in a statistically significant point source in an HET or 
LET image. This is the basis of our analysis. In this section we calculate the sensitivity for 
a point source in the center of the FOV of a single sub-telescope, and in §5 we consider how 
this sensitivity varies across the sub-telescope arrays' FOV. 

Skinner (2007) derived the source detection sensitivity of a coded mask system when 
standard assumptions are relaxed: the fraction of the open mask pixels may differ from 1/2; 
part of an open mask pixel may be occulted (e.g., by ribs around each pixel to support the 
EXIST masks); the detector pixels may not be small relative to the mask pixels; the source 
strength may be comparable to the background; and the closed mask pixels may be partially 
transparent (e.g., at high energy). Here we will consider the background-dominated case. 

Consider an image formed using counts accumulated over At and AE. The burst 
spectrum is N{E,t) (ph cm""^ s~^ keV~^). Let s be the source flux averaged over At and 
integrated over AE. Thus 




(1) 
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Let the average detector efficiency be 

_ lAEdEf^^dte{E)N{E,t) 
'° lAEdEl^,dtN{E,t) ' 

where e{E) the detector efficiency. At high energy photons leak through the closed mask 
pixels. Define an effective detector efficiency for fiux through the closed mask pixels 

_ lAEdEj^^dte{E)N{E,t)e-^-(^^ 

L^dEj^^dtN{E,t) ' 

where Tm{E) is the optical depth through the closed mask elements. Note that e includes 
absorption by all material over the entire detector, while Tjn{E) accounts only for absorption 
through the closed mask pixels. 

For the HETs Garson et al. (2006b,c) find that the photon aperture fiux from the Cos- 
mic X-ray Background (CXB) dominates the background below ~ 100 kcV, while at higher 
energy sources such as Earth albedo photons, charged particles and activation dominate the 
background. For the LETs the background results primarily from the CXB. The CXB con- 
tribution can be modelled semi-analytically, while other background sources require complex 
Monte Carlo calculations. Thus we model the total background count rate per detector area 
as 

B = (At)-' f dE [ dt [e{E)^E)Qa [/mask + (1 - /mask)e-"'"(^)] + b{E)] (4) 

Jae J At 

where /mask is the fraction of the mask area that is open, $(£') is the CXB (Cruber 1992), Via 
is the projected solid angle subtended by the detector's aperture averaged over the detector 
plane (calculated with the corrected formulae of Sullivan [1971]), and h{E) models the other 
sources of background. The aperture fiux includes the leakage of the CXB through the closed 
mask elements at high energy. 

The significance of the burst's image in the background-dominated case is 



S.^fJ^^^ . (5) 

where A is the detector area, and fm includes the factors resulting from the ribs around the 
mask pixels, the fraction of the mask pixels that are open, and the finite detector size. For 
open mask pixels where the ribs cover 0.2 of the pixel area and detector-to-mask pixel ratios 
of 1/2 (HET) and 1 (LET), /^ = 0.737 and 0.564, respectively 



It is convenient to parameterize the burst fiux in terms of 



Ft = (At)-' / dE ! dtN{E, t) = s ^ ° ^ ^ 
Jaeo J At Jae "-^ J At 



SAEjEj^ ,dtN{E,t) 
dtN(E, t) 



(6) 
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where A£;o=l-1000 keV. 

To understand the effect of transparency through the closed mask pixels, we define the 
mask factors 

1 _ e-'^'"(^) 
9s = 2 ' 

(relevant to the factor of (eo — ei)/2), and 

(relevant to the aperture flux). A decrease in Qs results in a decrease in the detection 
significance, while an increase in gi, indicates an increase in the aperture flux. 

The gamma-ray burst spectrum is modelled using the four-parameter 'Band' func- 
tion (Band et al. 1993): a low energy power law with an exponential roUoff {N{E) oc 
E°' exp[—E / Eq]) that merges smoothly with a high energy power law {N{E) oc E^). The 
break between the two power laws is characterized by Ep = {2 + a)EQ, which is the energy of 
the maximum of E'^N[E) oc vfy ii j3 < —2, i.e., Ep is the photon energy where most of the 
energy is radiated. We use the flux Ft (ph s~^ cm~^) integrated over the 1-1000 keV band to 
normalize the spectrum (see eq. 6). Thus the spectrum is characterized by the normalization 
Ft I the two spectral indices a and /3, and the energy Ep. 

Eq. 5 can be inverted to find the threshold value of Ft at the peak of the lightcurve 
for a given set of the spectral parameters that determine the shape of the burst spectrum — 
the spectral indices a and /3, and the peak energy Ep. The result is a surface in the four 
dimensional space given by these spectral parameters; bursts with spectra on one side of this 
surface (with Ft greater than the value on the surface) will be detected, while bursts on 
the other side will not. Holding the spectral indices a and (3 fixed projects this surface into 
a sensitivity curve in the Fx-Ep plane. The curve also depends on the accumulation time 
At (here At = 1). The dependence of the sensitivity on the accumulation time is discussed 
below (§6). 

Note that the threshold value of Ft at a given Ep is not the sensitivity of the detector 
at a photon energy equal to Ep. The power of sensitivity curves in the Ft-Ep plane is that 
they show the detectability of a burst with a given set of spectral parameters, and thus the 
sensitivity of different detectors can be compared, regardless of their specific energy response. 
A CZT-based detector that detects photons in the 10-150 keV band can be compared to a 
scintillator-based detector that detects photons in the 50-300 keV band. 
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4. Single Sub- Telescope Energy Sensitivity 
4.1. HET 

A single HET will have a 56 cm x 56 cm (an area of 3136 cm^) CZT detector plane that 
is 5 mm thick. The platinum and gold cathode pads on the CZT (~ 1000 A each), the Mylar 
thermal blankets (two ~ 5 mil blankets) and the Kevlar micrometeriod shield (~5 mil) in 
the current design produce negligible absorption > 10 keV. Figure 1 shows the efficiency of 
the CZT detectors as a function of energy 

We calculate the mask factors (eqs. 7 and 8) using the optical depth through a 5 mm 
thick, 107.9 cm x 107.9 cm plate of tungsten. Because of the supports necessary for the 
closed mask pixels, we assume an open fraction of /mask = 0.4. Figure 2 shows the resulting 
mask factors as a function of energy for both the nominal 5 mm thickness (solid curve) and 
for a mask with half this thickness (dashed curve). 

The linear dimensions of the detector and mask pixels will be 0.125 cm and 0.25 cm, 
respectively. For a detector to mask pixel dimension ratio of 1:2 the factor compensating for 
the finite size of the detector pixels is — 0.737 (see eq. 5). 

The background is modelled (eq. 4) as the sum of the CXB aperture flux and the 
continuum background from other sources (see Carson et al. 2006b, c for details). Figure 3 
shows the resulting background. 

We assume that counts are accumulated over two energy bands AE — 10-600 and 
40-600 keV. The 10-600 keV band is sensitive to soft bursts, where the large number of 
low energy burst counts compensates for the large aperture flux, while the 40-600 keV 
band is particularly sensitive to hard bursts where there are sufficient burst photons above 
40 keV. The required threshold significance is assumed to be 5*/ = 7 (the same as the BAT's 
threshold). 

Figure 4 shows the sensitivity curve for a single HET sub-telescope for three sets of 
spectral indices and for At=l s. 

The survey will localize sources to better than 56 arcsec (Crindlay 2007). The survey's 
threshold will be 5<t whereas the burst threshold will be 7a, and typically localization is 
proportional to ~ (cr — 1)~^. Thus we estimate that the HET's localizations will be better 
than 40 arcsec. 
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4.2. LET 

As currently designed, a single LET will have a 20 cm x 20 cm (an area of 400 cm^), 
1 mm thick Si detector plane with 0.02 cm pixels. Figure 1 also shows the LET efficiency. 
The mask will be 72 cm above the Si detector plane, with collimators extending from the 
detector plane to the mask. The 20 cm x 20 cm mask will have a thickness of 0.05 mm and 
0.02 cm pixels. Again, because of the need to support the closed mask pixels, we assume 
an open fraction of /mask=0.4. Over the LET's energy range (3 to 30 keV) the closed mask 
pixels will be optically thick. With a mask to detector pixel ratio of 1:1, the mask factor in 
eq. 5 is /^=0.564. 

Although we include an internal background of b{E) — 10~^ cts cm~^ keV~^ s~^ in 
our calculation, the background is almost entirely the result of the CXB aperture flux; see 
Figure 3. 

We assume a threshold image significance of 5"/ = 7 will be required over a single trigger 
energy band of AE— 3-30 keV. Figure 4 shows the resulting sensitivity. Note that the LETs 
are less sensitive than the HETs. The slopes of the HET and LET sensitivity curves are 
consistent with the energy dependence of the LET and HET detectors. 

The EXIST survey's locahzations should be better than 11 arcsec (Grindlay 2007), and 
thus accounting for the difference in survey and burst thresholds (5(7 vs. 7a), the LET's 
burst locahzations should be better than 8 arcsec. 



5. Off-Axis and Multi-Detector Sensitivity 

The arrays of HET and LET sub-telescopes will each cover a very large total FOV. 
Any point in these total FOVs will be in the fully- or partially-coded FOVs of a number of 
sub-telescopes. The resulting multi-detector sensitivity across the arrays' FOVs will depend 
on how the images from the different detectors will be added together; this merging will 
depend on the exigencies of the available computational power and the required data latency. 
Specifically, burst detection and localization on-board the EXIST spacecraft by radiation- 
hardened processors will probably be less sensitive than on the ground, where farms of 
high-speed processors will be available. In addition, localization on-board must be rapid so 
that telescopes on the ground can begin following the burst afterglow. 

The calculations above provide the on-axis sensitivity for single HET or LET sub- 
telescopes. Let R be the ratio of the actual sensitivity at a given point in the FOV to this 
single sub-telescope on-axis sensitivity, where sensitivity is proportional to Sj (see eq. 5) or 
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to the inverse of the threshold peak flux Ft- Thus larger R means a greater significance for 
a given peak flux or a smaller threshold peak flux for a given significance. 

The source flux faUing on the detector plane is only a fraction cos 9 of the flux it would 
have on-axis, where /c is the 'coding fraction' which accounts for the partial shadowing of the 
detector plane by the collimators (the detector sides) and 9 is the inchnation angle (the angle 
between the source direction and the detector normal). The non-source flux that contributes 
to the background around the source is proportional to the 'coding fraction' /c. In coded 
mask imaging only the counts in the region of the detector plane that is not shadowed for a 
given source contribute to the image around the source. The source flux that impinges on 
this region is foreshortened by the inclination angle (the 'cos 6'' effect), but the background 
in this region does not depend on the source's direction. For a single sub-telescope the ratio 
of the off-axis to on-axis significance is therefore R — fl^"^ cos 9. 

The methodology by which the data from multiple detectors will be combined is cur- 
rently being studied. First, the images can be added. Then the source flux is proportional to 
J2 fc,i cos 9i, the background to J2fc,i, and thus Rj = [Yl fc,i cos 9i]/ ^/YJ~. Alternatively, 
forming images for each sub-telescope and adding the significances for the common point 
sources in quadrature gives Rq — a/X] /c,i cos^ 9i. In practice for the HET and LET arrays 
the sensitivity over the FOV for these two methods differ very little, and we use Rj. 

To calculate the sensitivity over the FOV, we work, and plot results, in a coordinate 
system that is a projection of the spherical sky directly onto a plane peipeiidicular to the 
zenith, i.e., if a point on the sky has the coordinates x,y,z (where x'^ + y"^ + z'^ — 1), then 
we work in the x-y plane. In this coordinate system, z is along the spacecraft's zenith, x is 
along the direction of orbital motion, and the spacecraft nods (rocks) in the y direction. We 
calculate Rj at different points on this grid. 

Figure 5 shows the burst sensitivity over the sky for the HET and LET arrays; the 
maxima are just under twice the sensitivities (i.e., more sensitive than) of single HET and 
LET sub-telescopes. Figure 6 shows the amount of solid angle at a given sensitivity for both 
arrays. Thus different points in the overall FOV will have different sensitivity thresholds, 
which must be considered when analyzing the cumulative intensity distribution. Figure 7 
shows the low end of the cumulative intensity distribution resulting from variations in the 
threshold over the FOV; other effects that smooth the threshold are ignored, and therefore 
the effect demonstrated by this figure applies to any burst intensity distribution. Note that 
the sensitivity of Swiffs BAT also varies over the FOV, affecting the shape of the cumulative 
fluence or peak flux distributions. 
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6. Dependence on Accumulation Time At 

The HET and LET sensitivity curves presented in Figure 4 assumed At=l s, i.e., that 
the bursts were detected in images formed over 1 s. However, modern burst detectors 
(e.g., Swiff s BAT, the GLAST Burst Monitor and EXIST) usually use a number of different 
accumulation times. For an imaging detector the relevant At is the accumulation time for the 
final image. An accumulation time comparable to the burst duration will usually maximize 
the source significance. A longer accumulation time will dilute the signal with background, 
reducing the signal-to-noise ratio, and therefore the significance of the detection. On the 
other hand, a shorter accumulation time will often exclude signal that could have increased 
the significance of the burst detection. 

Quantitative analysis of the dependence of burst sensitivity on the accumulation time 
is difficult because of the large range of burst durations and the great diversity of burst 
lightcurves. Some bursts consist of contiguous, overlapping pulses while others have widely 
separated pulses. Band (2002) ran a software rate trigger with a wide range of At values on 
the lightcurves of 100 bright BATSE bursts, and determined that using a range of At values 
would increase the burst detection rate by ~25% over the rate for At=l s. Band (2006) 
explained the larger fraction of long duration bursts relative to short duration bursts in the 
Swift data set compared to BATSE's as resulting in part from Swiffs long accumulation 
times. 

As a demonstration of the increase in sensitivity afforded by using a variety of accumu- 
lation times, consider a burst lightcurve with an exponential shape, N{t) = A^o exp[— t/T]; 
the traditional duration of 90% of the emission is Tgo = T In 10. In this example the accu- 
mulation time is assumed to begin ait — O. Let Fr(At) be the threshold peak fiux averaged 
over 1 s (this is the quantity plotted in Figure 4 for At=l s) for a given At. Then the ratio 
of threshold peak fiuxes for two different accumulation times Ato and Ati is 

FrjAt,) ^ /At7 l-exp[-Ato/r] 

FriAto) V Ato 1 - exp[-Ati/r] ' ^ ^ 

If a detector uses a set of At values, then the smallest value of Ft (At) should be used for 
any given value of T. We assume we are in the background-dominated case (eq. 5); the 
detectability of very short bursts might be limited by a paucity of source counts. 

Figure 8 shows this ratio for Ato=l s and different sets of Ati. Thus this figure shows 
how the sensitivity of a mission such as EXIST to short and long duration bursts is increased 
by using a variety of accumulation times. The dashed curve assumes Ati=l s, and thus the 
ratio is equal to 1. Currently EXISTs planned imaging trigger (which is not the final imaging 
step in EXISTS burst detection process) will use Ati={3, 18, 108, 648 and 1296} s; this is 
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shown by the sohd curve. Finally, A^i may be varied to maximize the signal-to-noise ratio, 
minimizing Fr(Aii) to the smallest possible value; this is shown by the dot-dashed curve. 

If burst lightcurves could be described by the exponential shape of this example (and 
bursts did not undergo spectral evolution, which makes the duration energy-dependent), 
then the HET or LET threshold peak flux of a burst of a given peak energy Ep and duration 
Tqo would be the product of Ft from Figure 4 and the ratio from Figure 8. 

We emphasize that this is a highly idealized example meant to demonstrate how the 
variable accumulation times of EXISTs burst detection system will increase the sensitivity 
to long and short duration bursts. This is particularly relevant to high redshift bursts whose 
durations will be time-dilated. 



7. Discussion 

From the preceding analysis, we can draw several conclusions on EXISTs impact on 
the study of gamma-ray bursts. 

First we estimate the EXIST burst detection rate. The BATSE observations provide 
the cumulative burst rate as a function of the peak flux value ipB averaged over At = 1 s in 
the AE =50-300 keV band (Band 2002): 



Nh ~ 550 



/ 1 -0.8 



0.3 ph cm s 



bursts yr ^ sky ^ . (10) 



The HET threshold sensitivity for a single sub-telescope on-axis is ■0b ~ 0.12 ph cm"^ s"^ 
for Ep > 100 kcV. Using the BATSE rate in eq. 10 and integrating over the solid angle 
distribution in Figure 6 gives a burst detection rate for the HET of ~ 400 bursts per year. 
Note that this rate is over the BATSE-specific values of AE and At, and EXIST will use at 
least two different values of AE (see §4.1) and a variety of At values (see §6). Consequently 
this rate should be increased by approximately 50% to account for the soft, faint, long 
duration bursts to which BATSE was less sensitive than EXISTs HET will be; we therefore 
expect the HET array to detect ~ 600 bursts per year. 

The value oiipB for an LET varies more with the burst spectral parameters than for an 
HET, and therefore estimates of the LET burst detection rate based on the BATSE rate are 
much more uncertain. For a single LET iI^b ~ 0.3 ph cm~^ s"^ on axis at Ep — 100 keV, which 
gives a burst detection rate of ~180 bursts per year using eq. 10 and the LET distribution 
in Figure 6. This rate should be increased by a factor of 2 to account for the different energy 
band AE and accumulation times At. We use a larger adjustment factor for the LETs 
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than for the HETs because the LETs' energy band will overlap less with BATSE's than the 
HETs'. We therefore expect the LET array to detect ~ 350 bursts per year. 

Next we simulate the spectra that the EXIST suite of detectors will observe. Figure 9 
shows a count spectrum (counts s~^ keV~^) for a moderately strong burst as it might be 
observed by the LETs (lefthand set of curves), HETs (middle set) and the Csl active shields 
for the HETs (righthand set; based on Garson et al. 2006a). The solid curves show the 
signal count rate, while the dashed curves provide the estimated background. Thus EXIST 
will facilitate spectral-temporal studies. 

Particularly important to physical burst emission models is determining Bp, which is 
typically of order 250 keV (Kaneko et al. 2006). In addition, correlations of Ep with other 
burst properties, such as the 'isotropic' energy (the Amati relation — Amati 2006) or total 
energy (the Ghirlanda relation — Ghirlanda, Ghisellini & Lazzati 2004), have been proposed. 
'Pseudo-redshifts' calculated from the observables related to the burst-frame parameters in 
these relations can be used in burst studies when spectroscopic redshifts are not available, and 
can guide ground observers in allocating telescope time to observing potential high redshift 
bursts. The recently proposed Firmani relation (Firmani et al. 2006) correlates Ep, the peak 
luminosity, and a measure of the burst duration, all of which are related to observables in 
the gamma ray band. Thus pseudo-redshifts will be estimated using the Firmani relation 
based on EXIST data alone, independent of observations by other facilities. 

With well determined broadband spectra down to 3 keV, EXIST will be capable of 
determining whether the Band function (Band et al. 1993) suffices to describe burst spectra. 
For example, Preece et al. (1996) found evidence in the BATSE data for the presence of 
additional emission below 10 keV. 

By scahng from the EXIST survey's source localization (Grindlay 2007), we find that 
bursts should be localized at threshold by the HETs and LETs to better than 40 arcsec and 
8 arcsec, respectively; this localization should scale as ~ (c — 1)^^- Because the HETs are 
more sensitive to the LETs, the HET localization is relevant to the faintest bursts EXIST 
will detect. 

EXISTs burst capabilities calculated above will constitute a major leap beyond current 
detectors, and should increase the number of high redshift bursts detected. On average, high 
redshift bursts should be fainter, softer and longer than low redshift bursts (although the 
broad burst luminosity function and great variety in burst lightcurves and spectra obscure 
this trend). Figure 10 compares the detector sensitivities of the HET (solid curve) and LET 
(dashed curve) arrays to the BAT on Swift (dot-dashed curve) and BATSE's Large Area 
Detector (LAD — dot-dot-dashed curve). As discussed above, the sensitivity is the threshold 
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peak flux Ft integrated over the 1-1000 keV band as a function of the spectrum's Ep-, a — —1 
and f3 = —2 are assumed. In addition, the figure shows famihes of identical bursts at different 
redshifts (the curves with the points marked by '+')• Each family is defined by the value of 
Ep in the burst frame; here again a = —1 and f3 = —2 are assumed. In each family the burst 
would be observed to have Ft=7.5 ph cm~^ s~^ if it were at z = 1. The points marked by '+' 
are spaced every Az =1/2; thus the uppermost points are at z = 1 and the lowermost points 
are at z — 10. The pulses in burst lightcurves become narrower (shorter) at higher energy, 
an effect that is generally proportional to E^-^ (Fenimore et al. 1995). Since the observed 
lightcurve originated in a higher energy band, pulses should become narrower with redshift, 
reducing the peak flux when integrated over a fixed accumulation time; the plotted families 
include this effect. Finally, in §6 we showed that forming images on long timescales increases 
the sensitivity to long duration bursts, as might result from cosmological time dilation. 



8. Summary 

We presented our method for analyzing the gamma-ray burst sensitivity of EXIST, and 
applied it to the design for the Beyond Einstein program; this methodology will be used to 
guide and evaluate the evolving mission design. With two arrays of coded mask detectors 
covering the 3-30 keV (Si) and 10-600 keV (CZT) bands and non-imaging high energy Csl 
detectors (0.2-10 MeV), EXIST will be a significant gamma-ray burst observatory. EXIST 
will detect and localize ~ 2 bursts per day, observing their spectra from 3 keV to over 10 MeV. 
For bursts with comparable spectra and lightcurves EXIST will be approximately four times 
more sensitive than Swiffs BAT with a much larger FOV. With these capabihties, EXIST 
will accumulate a large sample of bursts with well determined properties such as Ep and 
redshift, facilitating physical modelling and population studies, and realizing the potential 
of gamma-ray bursts as cosmological probes. 
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Fig. 1. — Efficiency of the 5 mm thick CZT detectors (sohd curve) and the LET detectors 
(dashed curve) as a function of energy. 
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Fig. 2. — HET mask factors as a function of energy for a tungsten mask that is 5 mm thick 
(sohd curves) and 2.5 mm thick (dashed curves). The set of curves beginning at 0.4 at low 
energy is the mask factor gb (eq. 8), the fraction of the incident burst flux that reaches the 
detector plane, whether through an open or closed mask pixel. The set beginning at 0.5 at 
low energy is the mask factor gg (eq. 7), the fraction of burst flux that will be attributed to 
the burst when an image is formed. 
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Fig. 3. — The background spectrum for the HETs (sohd and dashed-dot curves) and LETs 
(dashed curve). The cosmic X-ray background (CXB) aperture flux (dashed curve for LET 
and dot-dashed curve for HET) is based on the parameterization of Gruber (1992) while the 
total HET background (sohd curve) adds other background components from Monte Carlo 
simulations (Garson et al. 2006b, c). 
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Fig. 4. — Maximum HET (lower set of curves) and LET (upper set) detection sensitivities 
for At = 1 s, the threshold peak flux over 1-1000 keV as a function of the spectrum's Ep. 
Solid hne — a — —1, P — —2; dashed hne — a — —0.5, P — —2; dot-dashed hne — a — — 1, 
P — —3. Note that this figure shows the sensitivity for detecting a burst with a spectrum 
characterized by Ep, and not for detecting a photon with an energy equal to Ep. 
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Fig. 5. — Burst detection sensitivity over the sky for the HET (left) and LET (right) arrays, 
in units of the sensitivity of a single sub-telescope on-axis. Greater sensitivity results in a 
smaller threshold flux. 



Fig. 6. — Solid angle at a given burst sensitivity for the HET (left) and LET (right) arrays, 
in units of a single sub-telescope. Greater sensitivity (larger value of R) results in a smaller 
threshold flux (smaller Ft). 
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Fig. 7. — Effect of the variable detection threshold across the FOV on the cumulative 
intensity distribution for the HET (solid curve) and LET (dashed) arrays. The assumed 
actual distribution (dot-dashed curve) is a power law with an index of -0.8. The intensity is 
given in units of the threshold value for a single sub-telescope. 
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Fig. 8. — Normalized burst sensitivity as a function of burst duration for different accumu- 
lation times At. The normalized burst sensitivity is the ratio of the peak flux averaged over 
1 s for At—1 s to the peak flux for different At values; a smaller value means the detector 
is more sensitive. The lightcurve is assumed to be an exponential in time. When a detector 
system employs more than one At then the minimum normalized sensitivity (resulting from 
the smallest peak flux) is used. The dot-dashed curve assumes At—1, and therefore is 1 for 
all durations. The solid curve shows the accumulation times currently planned for EXISTs 
image triggers: At=3, 18, 108, 648 and 1296 s. The dashed curve assumes all possible values 
of At and thus shows the smallest possible value of the normalized sensitivity. 
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Fig. 9. — Source count (solid curves) and background (dashed curves) spectra for the LETs 
(left hand set of curves), HETs (middle set) and Csl shields (right hand set). The burst 
spectrum has a. — —1, /3 = —2, — 300 keV and Ft — 7.5 ph cm~^ s~^. Based on 
the sub-telescopes' FOVs in the current design, we assume spectra can be formed from the 
equivalent of four HET and LET sub-telescopes, and the shields of nine HET sub-telescopes. 
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Fig. 10. — Maximum detector sensitivity for HET (solid curve), LET (dashed curve), Swiff s 
BAT (dot-dashed curve) and BATSE's LAD (dot-dot-dashed curve) assuming At = 1 s, 
a = —1, and f3 = —2. The HET and LET sensitivities assume burst detection by multiple 
sub-telescopes. Also shown are tracks for identical bursts at different rcdshifts. The bursts 
have different Ep and Ft=7.5 ph cm~^ s~^ when at z = 1. The points on the track arc spaced 
by Az =1/2; the faintest bursts on each track are at z = 10. Burst pulses are assumed to 
narrow by E"^-^. The assumed cosmology is = 0.3 and f^A = 0.7 {Hq is irrelevant to this 
calculation) . 
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Table 1. Parameters of the EXIST Detectors 



Parameter High Energy Telesco]:)e (HET) Low Energy Telesco]:)e (LET) 



Number 


19 


32 


Detector Material 


CZT 


Si 


Detector Thickness 


0.5 cm 


0.1 cm 


Detector Plane Dimensions 


56 cm X 56 cm 


20 cm X 20 cm 


Detector Pixels 


0.125 cm 


0.02 cm 


Mask Material 


Tungsten 


Tungsten 


Mask Thickness 


0.5 cm 


0.05 cm 


Mask Dimensions 


107.9 cm X 107.9 cm 


40 cm X 40 cm 


Mask Pixels 


0.25 cm 


0.02 cm 


Detector-Mask Distance 


140 cm 


72 cm 


Angular Resolution (FWHM) 


6.86' 


1.35' 


Localization {7a, 90% conf.) 


< 40" 


< 8" 


Fully Coded FOV 


21° X 21° 


16° X 16° 


Trigger Band AE 


10-600 keV 


3-30 keV 




50-600 keV 




/mask 


0.4 


0.4 


fm 


0.737 


0.564 



